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EXPERIMENTAL EVALUATION OF ANALYTICAL MODELS FOR 
THE INERTIAS AND NATURAL FREQUENCIES OF FUEL 
SLOSHING IN CIRCULAR CYLINDRICAL TANKS 
By Robert W. Warner and John T. Caldwell 


SUMMARY 


A correlation is presented of experimental measurement and analytical 
prediction of free, planar oscillations of a tank of mercury hung on a 
pendulum. Measured frequencies are presented for the slosh mode and the 
pendulum-tank mode, and time-history wave forms are included for qualita- 
tive comparisons. The sloshing Is represented analytically by a single 
pendulum or a spring-bob, and the size and location of their inertias are 
determined by matching forces and moments on the tank surface with line- 
arized hydrodynamic theory for inviscid fluids . It is concluded that such 
analytical models are identical for small displacements in the present 
application (as they are in simpler applications) and that they provide 
excellent simulation of the stiffness and inertia terms for fuel sloshing. 


INTRODUCTION 


A severe aero-servo-elastic stability problem can occur for vehicles 
which have a flexible structure and a swivelling rocket motor with high- 
gain control response. In an analytical study of such a problem fuel 
sloshing can be important as long as appreciable liquid fuel remains 
unburned. It is then necessary to determine whether the sloshing has a 
direct effect on stability, an indirect effect by modification of a defor- 
mation pattern, or no effect at all . Such an analytical study is often 
facilitated when the fuel sloshing Is incorporated in terms of analytical 
models, specifically pendulums or spring-bobs. Hence the main purpose of 
the present report is to evaluate experimentally the accuracy of these 
analytical models for a portion of the fuel sloshing problem, namely, the 
stiffness and inertia terms for sloshing in circular cylindrical tanks 
undergoing oscillations in a plane. A secondary purpose is to compare the 
pendulum and spring -bob models. 

A great deal of work has been done on sloshing of fluids in containers, 
and reference 1 is an excellent survey with an extensive bibliography. 
Guthrie (ref. 2) and Rayleigh (ref. 3 ) provided one of the earliest corre- 
lations between experiment and theory, showing excellent agreement as to 
fluid natural frequency. These results, however, are for stationary tanks. 
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An extensive correlation on the frequency response of the fluid force 
on the wall of an oscillating tank can be found in references 4 to 7 • 

Among these references there are comparisons of experiment and theory for 
high-damping tank configurations with baffles and low-damping configura- 
tions with bare tank walls . It is worthwhile to focus attention on the 
off-resonance correlations and the natural frequencies for the low-damping 
frequency responses (appearing in refs. 5-7) in order to evaluate errors 
in the stiffness and inertia simulation with minimum effect of the errors 
in the damping simulation. On that basis , references 5 to 7 indicate 
generally satisfactory stiffness and inertia simulation for the fluid 
force on the tank wall. Unfortunately, there are no moments measured on 
the surface of the tanks in references 5 to 7* 

The present experimental evaluation is based on low- amplitude , free, 
planar oscillations of a small, smooth-walled tank of mercury hung on a 
pendulum (the weight ratio of the mercury to the present rigid tank being 
representative of actual missiles, in which the fuel can be massive enough 
to have a significant effect) . This test setup is simulated mathematically 
by linearized equations of motion for an undamped two-degree-of-freedom 
system, with fuel sloshing represented by a single pendulum or a spring- 
bob. Solutions are compared with the experimental oscillations in terms 
of frequencies and wave forms. These measured and predicted quantities 
should be significantly affected by the fluid moments on the tank surface . 
Small amplitudes and linearized equations are used throughout for simplic- 
ity since oscillations universally begin at low amplitude unless an 
artificial initial condition is applied. 

Simple analytical models for fuel sloshing have also received a great 
deal of attention. Reference 8 shows that a system with spring -bobs and 
a fixed mass and rotary inertia exactly duplicates the tank-surface forces 
and moments given by linearized hydrodynamic theory for a rectangular tank 
oscillating in lateral translation, pitch, and yaw. References 6 and 9 
present similar exact systems for a cylindrical tank. 

The historical development of the pendulum model has not been as free 
of error as that of the spring-bob model. Reference 3 gives the correct 
lengths for pendulums which match the experimental and theoretical fluid 
natural frequencies. Reference 10 , however, presents pendulum masses 
based on duplicating the fluid force given by hydrodynamic theory on the 
wall of a rectangular tank restricted to zero motion. Unfortunately, the 
case of a stationary tank is not suitable for determining the pendulum 
masses, being useful only for finding the relative magnitudes of initial 
conditions for a fluid mode and its pendulum model after the masses have 
already been determined by another method. The correct pendulum masses 
and pivot-point locations for a cylindrical tank, plus the magnitude and 
location of a fixed mass and the magnitude of its rotary inertia, can be 
found among references 11 to 13 , as derived on the basis of forces and 
moments on the surface of a tank undergoing forced lateral and pitching 
oscillations . 
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As a result of the history just cited for the pendulum model, it 
appeared useful to compare the relationship of the pendulum model to the 
spring -bob model. For completeness, a brief development of the two models 
is given in the appendix prior to their comparison. 


NOTATION 


A infinite summation defined in equation (A5) 

a interior radius of tank 

b effective gravity 

eg center of gravity 


e 



■^xnpl 
^xns J 


-xop 

p 

p xos 



^6np\ 

■^0ns J 



Naperian base 

total fluid force on the wall of S. tank undergoing forced lateral 
oscillations, as determined by hydrodynamic theory and the 
pendulum model, respectively; positive right 

that part of the fluid force on the wall of a tank undergoing 
forced lateral oscillations which is contributed by the nth 
pendulum and the nth spring-bob, respectively; positive 
right, n = 1, 2, . . ., N 

that paid of the fluid force on the wall of a tank undergoing 
forced lateral oscillations which is contributed by the fixed 
mass in the pendulum model and the spring -bob model, 
respectively; positive right 

total fluid force on the wall of a tank undergoing forced 
pitching oscillations (about the center of gravity of the 
undisturbed fluid) , as determined by hydrodynamic theory and 
the pendulum model, respectively; positive right 

that part of the fluid force on the wall of a tank undergoing 
forced pitching oscillations (about the center of gravity of 
the undisturbed fluid) which is contributed by the nth 
pendulum and the nth spring-bob, respectively; positive 
right, n = 1, 2, . . ., N 

that part of the fluid force on the wall of a tank undergoing 
forced pitching oscillations (about the center of gravity of 
the undisturbed fluid'' which is contributed tu r the fixed mass 
in the pendulum model and the spring -bob model, respectively; 
positive right 



fluid height 


mass moment of inertia of the empty pendulum tank about its 
pivot point 

mass moment of inertia of the fixed mass about its own center 
of gravity in the pendulum model and the spring-bob model, 
respectively 

common symbol for I Q p and I os after they are proved equal 

total spring stiffness acting on the nth mass in the spring- 
bob model, n = 1, 2, . . . , N 

length of the nth simple pendulum in the pendulum model, 
n — 1, 2, • • •, N 

shortest distance from the hinge line to the center of gravity 
of the empty pendulum tank 

shortest distance from the hinge line of the pendulum tank to 
the center of gravity of the undisturbed fluid 

shortest distance, positive downward, from the hinge line of the 
pendulum tank to the hinge line of the pendulum in the single- 
pendulum model and the attachment line of the spring in the 
single-spring-bob model, respectively 

shortest distance, positive downward, from the hinge line of 
the pendulum tank to the center of gravity of the fixed mass 
in either the single pendulum or the single-spring-bob model 

shortest distance, positive downward, from the center of gravity 
of the undisturbed fluid to the hinge line of the nth 
pendulum in the pendulum model and the attachment line of the 
nth spring in the spring -bob model, respectively, 
n = 1, 2, . . ., N 

shortest distance, positive downward, from the center of gravity 
of the undisturbed fluid to the attachment line of the fixed 
mass in the pendulum model and the spring-bob model, 
respectively 

common symbol for Z 0 p and Z os after they are proved equal 

mass of the empty pendulum tank 

total fluid moment (about the center of gravity of the 

undisturbed fluid) on the surface of a tank undergoing forced 
lateral oscillations, as determined by hydrodynamic theory and 
the pendulum model, respectively; positive clockwise 
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^xnpl 

MxnsJ 


M xopl 

^XOSJ 



Menpl 

M0nsj 


M 0opl 
^0os J 


m 


that part of the fluid moment (about the center of gravity of 
the undisturbed fluid) on the surface of a tank undergoing 
forced lateral oscillations which is contributed by the 
nth pendulum and the nth spring -bob, respectively; positive 
clockwise, n = 1, 2, . . N 

that part of the fluid moment (about the center of gravity of 
the undisturbed fluid) on the surface of a tank undergoing 
forced lateral oscillations which is contributed by the fixed 
mass in the pendulum model and the spring -bob model, 
respectively; positive clockwise 

total fluid moment (about the center of gravity of the 

undisturbed fluid) on the surface of a tank undergoing forced 
pitching oscillations (about the same axis) as determined by 
hydrodynamic theory and the pendulum model, respectively; 
positive clockwise 

that part of the fluid moment (about the center of gravity of 
the undisturbed fluid) on the surface of a tank undergoing 
forced pitching oscillations (about the same axis) which is 
contributed by the nth pendulum and the nth spring -bob, 
respectively; positive clockwise, n = 1, 2, . . ., N 

that part of the fluid moment (about the center of gravity of 
the undisturbed fluid) on the surface of a tank undergoing 
forced pitching oscillations (about the same axis) which is 
contributed by the fixed mass with rotary inertia in the 
pendulum model and the spring-bob model, respectively; 
positive clockwise 

total mass of the fluid 


“npl 

E’nsJ 


“opl 

“osj 

mo 

N 


mass of the nth pendulum and the nth spring -bob, 
respectively, n = 1, 2, . . N 

common symbol for m^p and after they are proved equal, 

n = 1, 2, . . ., N 

mass of the fixed fluid mass in the pendulum model and the 
spring-bob model, respectively 

common symbol for mop and n^g after they are proved equal 

total number of pendulums or spring -bobs 

index indicating a specific mode or pendulum or spring -bob, 
n = 1, 2, . . ., N or « 


n 
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n g natural frequency of the slosh mode, that is, the mode which is 

mostly fuel motion, cps 

n-^ natural frequency of the pendulum-tank mode, that is, the mode 

which is mostly tank motion, cps 

s n lateral displacement of the mass of the nth spring "boh 

relative to the tank axis, positive right, n = 1, 2, • . ., N 

s xn the s n type motion which results from lateral oscillations 

of the tank, n = 1, 2, . . ., N 

the s n type motion which results from pitching oscillations 
of the tank about the center of gravity of the undisturbed 
fluid, n = 1, 2, . . . , N 

t time, sec 

x lateral displacement of the tank in lateral translation; 

positive right 

Xq amplitude of x in harmonic oscillation 

e n nth zero of ^‘(e), where Ji(e) is the first-order Bessel 

function of the first kind and the prime indicates differen- 
tiation with respect to the argument, n = 1, 2, . . ., oo 

9 pitching displacement of the pendulum tank about its hinge line, 

or of the tank alone about the center of gravity of the 
undisturbed fluid 

0 Q amplitude of 6 in harmonic oscillation 

pitching displacement of the nth pendulum relative to the 
tank axis, positive clockwise, n = 1, 2, . . ., N 

cp^ the Cf^ type motion which results from lateral oscillations of 

the tank, n = 1, 2, . . . , N 

cp^ the Cf^ type motion which results from pitching oscillations 

of the tank about the center of gravity of the undisturbed 
fluid, n = 1, 2, . . ., N 

w circular frequency of forced harmonic oscillation, radians/sec 

w n natural circular frequency of the nth fluid mode, the nth 

pendulum, or the nth spring -bob, radians/sec, n = 1, 2, 

. . . , N or oo 
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EXPERIMENTAL AND ANALYTICAL METHODS 
Apparatus and Test Technique 


The experimental data were obtained with the model shown in figure 1. 
Basically it consisted of a 3 .89-inch-inside-diameter tank mounted in a 
rigid frame and pivoted about a horizontal axis. The support arms could 
be moved along the longitudinal tank axis to vary the distance from the 
pivot point to the center of gravity of the empty tank. Mercury was used 
as the sloshing liquid. 

Fluid surface displacement at the tank wall was measured with a 
cylindrical capacitive probe partially submerged in the liquid, and with 
an accompanying exiernal bridge circuit. This measurement gave a time 
history of the fluid motion relative to the tank motion. Data for the 
tank motion were obtained with a noncontacting optical system. In this 
system a light beam is focused at a fixed point on the tank and through 
a servo control remains locked on the point for small tank displacements. 
The displacement of the light beam, which was proportional to the tank 
displacement, then activated a photoelectric cell and the output from the 
latter was recorded on an oscillograph. Similarly, the capacitive bridge 
circuit output, which was proportioned to the relative fluid motion, was 
recorded on an oscillograph. 

In general, oscillations were initiated by pulling the tank to one 
side, allowing all fluid motion to cease, and then releasing the tank. 

For accuracy in determining the slosh mode frequency, the mode was excited 
by manually forcing the tank at a frequency near that of the slosh mode. 
After excitation had been achieved, forcing was stopped and the subsequent 
free oscillation motion of the fluid was recorded. Under these conditions 
the fluid motion was almost entirely that of the slosh mode and thereby 
could be used to determine the slosh mode frequency. To eliminate system 
noise error for this case the bridge output was subjected to a frequency 
filter before it was recorded. 

For both the slosh mode and the pendulum-tank mode, the frequency for 
each pivot-point location and fluid height was determined three times. 
Repeatability was within 1 percent or better. Since an average of the 
three readings was used, the experimental error is estimated to be l/2 
percent or better. 


Range of Data and Sources of Error 


To check the effect of fluid height on experimental-analytical 
correlation, tests were performed over an h/a (ratio of fluid height 
to tank radius) range of 0*50 to 2.00. For each h/a value an l 
(distance from pivot point to empty-tank eg) range of roughly 0.6 to 2.5 
inches was tested to indicate the relative effects of rotation and 
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translation upon the correlation. For the smallest Z and largest h/a 
values tested (Z = 0.59^ in*; h/a = 2 .00) , a l/lO-inch transverse dis- 
placement of the undisturbed fluid center of gravity corresponded to an 
angular rotation of 3-77°* For the largest Z and smallest h/a value 
tested (Z = 2.469 in . , h/a = 0-5Q) the same transverse displacement 
corresponded to an angular deviation of 0*90°. Free oscillations were 
started from an initial pendulum-tank angle of around 10 , which was small 
enough to prevent splashing . 

The possible sources of error pertain, of course, to the quantities 
used to correlate experimental and analytical results, namely, the fre- 
quencies of the slosh mode (mostly fuel motion) and the pendulum-tank mode 
(mostly tank motion) , together with the wave form of the two measured 
coordinates . 

One source of error arises from the fact that damping was neglected 
in the analysis . Hence only qualitative wave form correlation would be 
expected on the time histories. On the other hand, in analogy to a linear 
one-degree-of-freedom system, slight damping would not change the modal 
natural frequencies appreciably. A computation using the observed damping 
factor indicated a natural frequency change of less than 1 part in 10 4 , 
well within the estimated maximum experimental error of 1/2 percent. 

The capacitive probe contributed two possible errors in the 
determination of the fluid motion, namely noise and the tendency of the 
fluid to stick to the probe. These effects increased the difficulty 
experienced in wave-form correlation for the fluid motion but had no 
effect on the frequency determination. 

The possible frequency error from recorder paper speed was checked 
by a precision frequency generator. Indicated error from this source was 
less than 1 part in 2000. This check was made for both the analog 
computer and experimental data recorders. Possible frequency error from 
filtered traces was determined by comparison with nonfiltered traces where 
possible. All such checks indicated deviations within the estimated maxi- 
mum experimental error. In a similar fashion manually excited slosh mode 
frequencies were checked against the same frequencies as determined from 
quiescent initial conditions . Again deviations were within the outer 
limit of 1/2 percent arrived at from the repeatability check and the 
averaging procedure . 


Analytical Models and Methods 


In the appendix it is shown that both the pendulum model and the 
spring-bob model can be made to duplicate exactly the tank-surface forces 
and moments due to forced lateral translation and pitching of a partially 
filled cylindrical tank (shown in fig. 2) as predicted by linearized 
hydrodynamic theory. These models (pictured in figs. 3(a) and 3(b)) can, 


9 


of course, "be used to determine parameters for the analytical simulation 
of free vibrations of the pendulum tank considered in the present report. 
Figure 4(a) shows this tank simulation for the case where a single 
pendulum of length Lx and mass m 1; and a fixed mass dIq with rotary 
inertia I 0 , are used as a fuel model. The remainder of figure 4(a) is 
self-explanatory except for the fact that the figure utilizes the following 
identities between the pendulum and spring-bob models taken from equations 
(A19) , (A20) , (A22) , and (A23) of the appendix for n = N = 1: 

m is J m ip ® m i 4os = ^op = ^o 

Hbs = “op = ™o Iqs = lop = Io 

where the symbols are defined under Notation. It should be noted that 
Cfh is proportional to the relative fluid motion in the first mode . 

With b for effective gravity, and M and I denoting the mass and 
hinge-line moment of inertia of the pendulum tank without fuel, the 
following two linearized equations of motion can be derived for the system 
of figure 4(a) by Lagrange's equations or direct application of force 
balance: 


L 1 CP 1 + ( + Lx) 0 + b(<Pi + 0) = 0 (l) 

(i + I Q + hIq 1-^ + m. x lf^ s )0 + m 1 ZppLi(cPi + 0) + b(MZ + Z^ + m x z fp )e - 0 

(2) 

Figure 4(b) shows the case where a single spring-bob of stiffness 
k x and mass m x replaces the pendulum in the fuel model. The two 
linearized equations of motion are as follows: 

* s i _ ®^fs + — '■ Sx - "b© = 0 (3) 


(i + I 0 + n^Z-t 2 + mxLpg 2 )© - mi2^ g sx + b(MZ + n^Z^ + mxZ.p s )0 - mxbsx = 0 




where sx is the relative displacement shown in figure 4(b) . 

Since the pendulum and spring-bob models must produce the same tank 
surface forces and moments on a mode-by-mode basis in order to duplicate 
the same hydrodynamic theory, it is expected that the linearized equations 
(3) and (4) can be converted to equations (l) and (2) by a simple linear 
transformation between Sx and cpi* This expectation is borne out when 
the substitution 


Sx - -LxCPi 


(5) 
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is made in equations (3) and (4) , together with the identities ( API ) , 

(a8) , (Al4) , and (A19) from the appendix for n = N = 1. Hence equations 
(l) and (2) were used interchangeably with equations ( 3 ) and (4) for the 
analysis . Occasionally both sets were used as a check on the analytical 
frequencies and wave forms. 

After the parameters for equations (l) and (2) or equations ( 3 ) and 
(4) were specified numerically as described in the next subsection, 
analytical frequencies were determined by extracting roots from the 
characteristic determinant of these equations. For convenience in plotting 
the time histories, the equations were also set up on an analog computing 
machine and subjected to an initial displacement in Q, with a correspond- 
ing initial Cfh. or si determined by static equilibrium. These initial 
displacements, together with zero initial velocities, correspond to the 
experimental initial conditions described in a previous subsection. 


Numerical Determination of Analytical Parameters 


The numerical parameters of table I were used for equations (l) and 
(2) (and various identities from the appendix for n = N = 1 would give 
immediate conversion if parameters for equations ( 3 ) and (4) were desired) . 

The value given in table I for the fixed mass M of the empty 
pendulum tank was measured directly, of course, and b was taken as the 
sea-level acceleration due to gravity. Different fluid heights h and 
various pivot-point locations were selected to determine the two 
parameters, h/a and Z, which are considered independent variables in table 
I. The tank radius a was measured directly, and Z was determined by 
balancing the empty pendulum tank horizontally on a knife *edge for each 
pivot-point location and measuring the distance from the pivot point to- 
the knife edge. For each value of h/a, the corresponding value of Li 
was determined from equation (a 8) in the appendix with e! = 1.84119 
(ref. l4) . Also for each h/a, the total fluid mass m was measured, 
and this permitted the determination of m x , m,-,, and I G by equations 
(A9), (AT), and (A13), respectively. For each value of Z, the hinge-line 
moment of inertia I of the empty pendulum tank was determined by its 
measured natural frequency. Finally, for each combination of Z and h/a, 
the distance Zj^ from the hinge line of the pendulum tank to the center 
of gravity of the undisturbed fluid was measured, and this was used with 
equations (AlO) and (All) and the length identities indicated in figure 
4(a) to determine the values of Z^ and Z^ given in table I. 
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RESULTS AML DISCUSSION 
Comparison, of Frequencies 


A total of 28 combinations of Z and h/a were used in this 
investigation. The slosh -mode and pendulum-tank-mode experimental and 
analytical frequencies are presented in table II. For this set the maxi- 
mum deviation of an analytical from experimental modal frequency was 2.9 
percent. Ninety percent of the frequencies had deviations of 1 percent 
or less; 70 percent had deviations of 0.5 percent or less; and percent 
had deviations of 0.2 percent or less. 

It should be noted that the frequencies discussed above are for the 
normal (combined) modes of a two -degree -of -freedom system in which coupling 
between the pristine modes (as designated herein) occurs. For the present 
system the latter consists of one mode in which the tank is held fixed 
while the fluid is in motion and a second mode in which the tank moves 
with the fluid frozen inside. What has been called the slosh mode is 
actually a combination of the two pristine modes in which the moving fluid 
component is predominant. Similarly, the pendulum-tank mode is a combi- 
nation in which the tank plus frozen fluid component predominates. 

Since the pristine modes are independent of all the fuel model 
parameters except those giving fluid natural frequencies, any free oscil- 
lation test of all the parameters requires that the coupling (resulting 
from fluid-tank force and moment interactions) described above be present. 
Furthermore for the test to be adequate the coupling effects must be 
sizable. The simplest measure of coupling in a system is the change in 
coupled modal frequency from pristine modal frequency. Computations 
revealed that for the present test such frequency changes ranged from 0.5 
to 4 !+ percent, with a roughly linear distribution for intermediate states. 
This indicates that enough data points with sizable coupling effects have 
been included to make the test significant; and the test then has the 
special advantage of determining how the analytical sloshing parameters 
perform when a combination of modes is present, as in an actual missile. 


Comparison of Time-History Wave Forms 


For a qualitative comparison of time-history wave forms, three sets 
of unfiltered analytical-experimental traces are shown in figure 5* These 
are representative of the various qualitative features used in the corre- 
lation. In figure 5(a) the 0 trace exhibits beating, with experimental 
and analytical beat frequencies nearly equal. The corresponding traces 
for cp 3 _ and fluid Ins igivt cLo n ot exhibit features as outstanding as the 
smooth beating just described. The wave-form correlation here is poor- 
It should be noted that the experimental traces for fluid height required 
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larger amplitudes of tank motion than did the corresponding • 9 traces. 

It is possible that these amplitudes introduced nonlinearities which, with 
the damping of the system and the limitations of the capacitive probe, 
made wave form correlation poor or indecisive for all fluid traces except 
the most outstanding. 

An example of outstanding wave-form correlation between s* (note in 
eq. ( 5 ) that S 3 . = -Li9i) and fluid height is shown in figure 5(b) . Here 
the easily correlated feature of one frequency component visibly super- 
posed on the second is present . The result is a series of irregular peaks 
and valleys that can be followed on both the analytical and exper ime nt al 
traces for several cycles. The observable noise has been discussed 
earlier. The accompanying 9 coordinate is featureless in this instance. 

Figure 5(c) represents a state of the system intermediate between the 
preceding two. Here the 9 trace exhibits "snaking," or more precisely, 
a low-frequency oscillation of the mean line for the most immediately 
apparent sinusoidal wave. The wave form agreement here is quite good. 

The corresponding Cff and fluid height traces are not particularly 
distinguishable, but both do display the same general peak-to-peak 
variations . 

On the average, an analytical time history duplicated the 
corresponding experimental coordinate time history for 4 cycles, with 
general agreement continuing indefinitely. 


CONCLUSIONS 


An experimental evaluation has been presented for specific 
analytical models for the stiffness and inertia terms of fuel sloshing in 
circular cylindrical tanks . The evaluation is based on two-degree-of- 
freedom analytical predictions, with sloshing represented by a single 
pendulum or spring-bob, of the experimental frequencies and wave forms for 
small, free oscillations of a smooth-walled pendulum tank containing 
mercury. The following conclusions are indicated: 

1- In view of the maximum error between experiment and analysis of 
2.9 percent for all slosh and pendulum tank frequencies, the maximum error 
of 0.5 percent for 70 percent of those frequencies, and the surprisingly 
good correlation of the time-history wave forms despite the damping of the 
physical system, it is concluded that either the pendulum or the spring-hob 
model provides excellent simulation of the stiffness and inertia terms for 
fuel sloshing . 
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2 . The well-known identity between pendulums and spring-bobs for 
small oscillations is maintained in their rather complicated simulation 
of tank-surface forces and moments given by linearized hydrodynamic 
theory for a partially filled cylindrical tank undergoing small lateral 
and pitching oscillations , and this simulation is exact. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 1 6, 1961 
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APPENDIX A 

DEVELOPMENT AND COMPARISON OF ANALYTICAL MODELS FOR 
FORCED LATERAL AND PITCHING OSCILLATIONS OF A 
PARTIALLY FILLED CYLINDRICAL TANK 

FORCES AND MOMENTS ON TANK SURFACE 


The method for developing analytical models for forced lateral 
translation and pitching of a partially filled cylindrical fuel tank is 
to determine masses and locations of pendulums or spring-hohs so as to 
match forces and moments on the tank surface, which are known from exact 
linearized hydrodynamic theory for inviscid fluids. The exact results, 
which can he derived by application of linearized boundary conditions to 
the Laplace equation as in reference l4, are presented below, with symbols 
defined under Notation and in figure 2: 



(Al) 


(A2) 


(A3) 
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M 0 = mat) 


4 * 

a 


0 o e iu)t + u> 2 0 o e ia)t ma 2 




e n coBh (e 




(AM 


where 


A 

A h 


+ 5 coshfe n - 


\ 


a n=i % 3 (% 2 - l)slnh (e n | 


(A5) 


and the pristine slosh frequencies u) n are given by 


“n 2 = I e n tanh (e n |) , n = 1, 2, 


\ 


(a6) 


These particular arrangements of the exact hydrodynamic results are 
important for matching with pendulums or spring -hots . For that purpose 
to should enter within a summation only in modal coefficients of 
l/[(w 2 /cj 2 ) - l], and everything which is not such a coefficient should 
be taken out of the summation and regarded as a summed number . Then the 
modal coefficients can be compared with those of the pendulums or spring- 
bobs on a mode-by-mode basis. 


FENDUIUM MODEL 


The pendulum model is shown in figure 3(a). For a total of N 
pendulums, giving N fluid modes, the following requirement is imposed: 


m = mop 



(A7) 


n=i 
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Deter mi nation of the model requires solution, for the unknowns 1^, n^p, 
mop , Z n p , 2 0 p , and I Q p • With the well-known formula for the natural 
frequency of a simple pendulum, equation (a 6) yields immediately 

= ± , n = 1, 2, . . N (A8) 

e n tanh fe n 



The first major step in solving for the remaining unknowns is to 
match F x an d M x for forced lateral motion of the tank of figure 3(a) 
at x = x^e^ 10 ^, positive right. The linearized equation of motion for 
m^p is (with cp^ defined in the Notation) 

%ip ( L n* c Pxn + ^xoeiwt) + nfcpbq^ = 0 

and solution for cp gives F^p, the contribution of the nth pendulum 
to the wall force p , by 

F xnp = -%p b< Pxn> n = 1, 2, . . N 
It is evident that for the fixed mass the force on the tank wall is 


F xop “op 1 ^ 2 ^ 6 


iwt 


and after a little algebra is given by 

N 


F xp “ F xop + / i F xnp 
n=i 



w. 


h 


- 1 


A mode-by-mode comparison of with F x of equation (Al) gives 


2m tanh (s n — 


“np 


e n ( € n 2 - !) 7 


(A9) 


a 


and equation (A7) can be used to solve for m^ . 

It is apparent from figure 3(a) that for forced translation the 
tank-surface moment Mjjp must have the following components: 

^xnp = _ ^np F xnp 

^xop = ~^op F xop 


Then 
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N 


Mjq? - M XO p +^M xn p 
EL— 1 

and when the algebra is worked out,, a mode-by-mode comparison of M™, 
with M x of equation (A2) yields 


l 


BP = 

h 


1 + 

2 


% | sinh(e n 1 


-n 


a 


tanh(€ n — 


(A10) 


Comparison of the l/|>(h/a)] term of equation (A2) with the correspond i ng 
term of yields 


l, 


°P = 
h 


m 


^np 

“op B mop 


n=i 


a 


(All) 


It is of some interest to note the sense in which equation (All) 
preserves the center of gravity of the fluid. If from equation (a8) 

is put in ) (mnp/mop) (L n /h) and this quantity is added and subtracted 
n=i 

on the right-hand side of equation (All ) , some algebraic manipulation 
yields 


N 


“op^op + ) m np(^np + In) = Iaa - 
n=i 


1 


* / e n 2 K 2 " 1) 


If N -*• it can be shown that ) = and thus in figure 

n=i £ i 2 ( £ i 2 - !) ° 

3(a) the pendulums preserve the center of gravity of the fluid. Even for 
a single pendulum (N = l) the center of gravity is virtually maintained 
since the value of e n for n = N = 1 yields l/[ei 2 (ej. 2 - l) ] = l/8.1. 


The second and final major step in developing the pendulum model of 
figure 3(a) is to match F q and M 0 for forced pitching of_the tank about 
the center of gravity of the undisturbed fluid at 6 = positive 

ClOGiCVTiSC • j7l1S i_XilG3,iTiZSu. 0 C^uS/t) 2. Oil. Ci luO’ti.Gii ±OT J-fc* \wxl-ii 

defined in Notation) 
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%pt L n'P©n “ w2 ( ^np +Ln)e o e iwt ] + 'bi% 1 (cp eri + ©oe^t) = 0 (A12) 

Again solution for cpfo gives Fg^, the contribution of the nth pendulum, 
to the wall force Fgp, by 

F ©np = — • m np"k < -P©n’ n = Ij 2, . . • , IT 

A combination of d’Alembert and gravity forces gives the force on the 
tank wall for the fixed mass as- 

F ©op “ m op l30 o elWt - m op z op w20 o eia,t 


F 0p F 0op ^ ^©np 
n=i 

and substitution of equations (A9) , (A10) , and (All) gives an exact check 
of Fg in equation (A3) for all terms outside the modal summation and 
an exact mode-by-mode check within the summation. This result is strictly 
a check, of course, since no new unknowns are involved. 


Figure 3(a) indicates that for forced pitching the tank-surface 
moment Mgp must have the following components: 


M ©np = - ^np F 0np 


M eo P " - Veop + “% P e o el “ t 


The use of equations (A9) , (lO) , and (All) in 

IT 

M©p — Mg 0 p + ^ 


n=i 


then gives an exact mode -by-mode check within the modal summation of 
Mg, equation (A4), and also checks the term mab jl/[4(h/a) l}©,^^ 0 ^. The 
remaining terms are matched by the following solution for lop: 

IT 

^op ~ + ^ “op^op 2 ~ ^“np^np 2 + (A13) 

n=i 

where A is the infinite summation of equation (A5) • It should be noted 
that the first two terms of I Q p give the moment of inertia of the 
fluid as if it were frozen. The solution for lop? in terms of the 
quantities i%p, H^p, ^np? and 4op found earlier, completes the determi- 
nation of the pendulum model of figure 3(a) . These results agree with 
previous results found among references 11 to 13 when the latter are 
extended beyond N = 1. 
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SPRING -BOB MODEL 


The spring-bob model is shown in figure 3(b), and equation (A7) with 
subscripts s instead of p, is imposed as for the pendulum model. 
Determination of the spring-bob model requires solution for the unknowns 
k n , ninsj mos> ^ns> Z os , and Ios* When the formula for the natural fre- 
quency of a spring -bob is combined with equation (A 6) , k n can be found 
immediately in terms of mng as follows 


k n = %is | € n tanh (e n n = 1, 2, 


N 


(AllO 


The procedure in solving for the remaining unknowns is identical to 
that for the pendulum model. The following key equations in the develop- 
ment of the spring model replace corresponding equations in the pendulum 
development (see Notation and fig. 3(b) for definitions of symbols): 

®ns( s xn - w XQel^l) + kj^s^^ = 0 

F xns = k n s xn 

Fxos = mos w2x o elu)t 


^xns = ^ns F xns + “ns^xn 
^xos _ ~^os F xos 

“ns( ¥ 0n + ^InsV 1 ^) + Vfci “ “ns 150 


0ns 


- V 


0n 


F 0OS = m os b0 o elWt - m os z os u20 o elwt 

M 0ns = " 2 ns F 0ns + m ns bs 0n 

M„ = -2 nc ,F flAt , + w 2 I 0 
0OS OS 0OS OS 0 


From these equations the procedure followed for the pendulum model yields 
the following results to supplement equation (A7), with subscripts s 
instead of p, and equation (Alb) for the spring -bob model: 


“ns 


2m tanh 




h 


Sn^n 2 - x ) a 


/,,r\ 
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(Al6) 


(ait) 



As in the case of the pendulum model, the foregoing results for the 
spring-hob model duplicate the exact hydrodynamic forces and moments on 
the tank surface even for those terms whose duplication requires the use 
only of previously determined unknowns . When N -> the spring-bob 
results are identical to those of references 6 and 9* By the argument 
immediately following equation (All), it is apparent that equation (A17) 
automatically preserves the fluid center of gravity as N -* °° and very 
nearly does so for a single spring (N = l) • 


RELATIONS BETWEEN PENDULUM AND SPRING-BOB MODELS 


Comparison of equations (A9) and (A15) shows 

®ns = ®np = “n > ^ — j ^ (A19) 

where 1 % is the symbol to be used for both models. Hence, by equation 
(A7) with subscripts p or s, 


“os = “op = “o 


(A20) 
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Comparison of equations (A10) and (Al6), and use of equation (a 8) , yields 

^ns = ^np + (A2l) 

When equations (A2l) , (A8) , and (A15) are introduced in equation (A17) , 
it is found that 

^os ~ ^op s ^0 (A22) 

where Z 0 is the new symbol for both models. A similar operation on 
equation (Al8), which involves lengthier algebra, yields 

■^os ~ ^op = ■*■0 (A23) 

None of the above exact comparisons requires N -* °o. In fact the 
n = 1 term of all summations not involving the modal coefficient 
l/[(m n 2 /w 2 ) - 1] is completely dominant except as h/a -*■ 0. As h/a -*• 0, 
the spring-bob forms in equations (Al6) through (Al8) are better for 
computational purposes than the corresponding pendulum forms since the 
hyperbolic sines and tangents in the denominator combine to give a half 
angle hyperbolic tangent in the numerator. For equation (Al8) this 
combination requires N -* °° since A of equation (A5) is an infinite 
summation. These computational matters are not too vital, however, since 
as h/a 0 the fuel sloshing effect on stability disappears. 
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TABLE I.- ANALYTICAL PARAMETERS 
[Symbols are defined in the section on Notation] 



h/a 

Li, 

in. 

m i» 

lb sec 2 /in. 

lb sec 2 / in. 

Io, 

lb-in. sec 2 

l/2 

1 

3/2 

2 

1.454 

1.110 

1.064 

1.058 

0.00971 

.01271 

.01325 

.01334 

0.00499 

.01668 

.03085 

.04542 

0.00684 

.00654 

.00952 

.02202 


l, 

in. 


0-594 

.969 

1.219 

1 .469 
1.844 
2.219 

2.469 


I, 

lb-in. sec 2 


0.1174 

.1264 

.1256 

.1285 

.1473 

.1668 

.1813 


\ h/ a 

1/2 

1 

3/2 

2 

Z \ 


z t> 

l 

h> 

Z fp> 

V 

4fp> 


in .\ 

in . 

in . 

in . 

in . 

in • 

in . 

in • 

in . 

0.594 

2.448 

2 .664 

2.447 

2.569 

1.849 

2-337 

1.027 

2.006 

.969 

2.948 

3.164 

2.947 

3.069 

2.349 

2.837 

1.527 

2.506 

1.219 

3.448 

3.664 

3.447 

3.569 

2.849 

3-337 

2.027 

3.006 

1.469 

3-948 

4.164 

3.947 

4.069 

3-349 

3.837 

2.527 

3.506 

1.844 

4.448 

4.664 

4.447 

4.569 

3.849 

4-337 

3-027 

4.006 

2.219 

4.948 

5-164 

4.947 

5.069 

4.349 

4.837 

3.527 

4.506 

2.469 

. 

5.448 

5 .664 

5.447 

5.569 

4.849 

5-337 

- . 1 ■ J 

4.027 

5.006 


































- NATURAL FREQUENCIES FOR NORMAL (COMBINED) MODES 
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Figure 1 . - Apparatus . 










(b) Spring-bob model. 

Figure 3 *- Concluded. 






31 











Ordinates are fluid height change at tank wall relative 
to tank bottom for the experiment and for 

the analysis. 




Ordinates are 9 for experiment and analysis. 

(a) h/a = I.OO, l = 2.219 inches 

Figure 5*- Sample experimental and analytical wave forms. 






Ordinates are fluid height change ot tank wall relative 
to tank bottom for the experiment and for 

the analysis. 




Ordinates are 9 for experiment and analysis 

(c) h/ a = 1.50, l = 0.969 inch 


Figure 5.- Concluded. 
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